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Abstract. Supersymmetric quantum mechanics is well known to 
provide, together with the so-called shape invariance condition, an 
elegant method to solve the eigenvalue problem of some one-dimensional 
potentials by simple algebraic manipulations. In the present paper, 
this method is used in statistical physics. We consider the local 
critical behaviour of inhomogeneous Ising models and determine the 
complete set of anomalous dimensions from the spectrum of the 
corresponding transfer matrix in the strip geometry. For smoothly 
varying perturbations, the eigenvalue problem of the transfer matrix 
indeed takes the form of a Schrodinger equation, and the corresponding 
potential furthermore exhibits the shape invariance property for some 
known extended defects. In these cases, the complete spectrum is derived 
by the methods of supersymmetric quantum mechanics. 
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1. Introduction 

The concept of supersymmetry has first appeared in quantum field theory and has 
been later used in different areas of physics (cf random systems) . The essence of the 
method is more transparently seen in ordinary quantum mechanics as it is known 
since the work of Witten [1]. Supersymmetric quantum mechanics [SSQM] provides a 
unified framework to perform the factorization of the Schrodinger equation, following 
the pionneering works of Dirac and of Schrodinger ([2-5], for a review see reference [6]). 
Furthermore, if the potential in the Schrodinger equation has the property of shape 
invariance [7] , the eigenvalues and the corresponding eigenvectors can be obtained by 
simple algebraic manipulations and it was found that the well known exactly solved 
problems (i.e. those problems which can be rewritten as hypergeometric equations 
after a suitable change of variable) exhibit the shape invariance property. 

In the present paper, we show a possible new field of application of SSQM. In 
statistical physics, inhomogeneous systems have been extensively studied in the past 
decade (for a recent review, see reference [8]). An inhomogeneity can be caused 
basically in two different ways. Geometrical effects due to the surface shape of the 
system and/or modified couplings or defects may influence the critical behaviour. The 
simplest inhomogeneity is the semi-infinite system with a free surface. The universal 
behaviour in a surface layer with a width of the order of the correlation length is 
described by a set of local (surface) critical exponents which are different from the 
bulk ones (see in reference [8]). More generally, the existence of a free surface may 
induce a coupling enhancement between nearest neighbour spins in a region of some 
extent close to the boundary and a local modification of the critical behaviour can 
then be expected. One special type of extended defect was introduced by Hilhorst 
and van Leeuwen [9]. Here the couplings perpendicular to the surface deviate from 
the bulk one as a power law A/y u , y being the distance from the free surface. It 
follows from a relevance-irrelevance criterion [10,11] that this type of perturbation 
is marginal for the two-dimensional Ising model at ui = 1. In this case the critical 
exponents are A-dependent as obtained from a number of exact calculations [12-18]. 
The conformal properties of such systems have been investigated using the plane-to- 
cylinder conformal mapping, under which the system is mapped onto a strip. Provided 
the perturbation profile is also properly transformed, the gap-exponent relation [19] 
and the tower-like structure of the spectrum are preserved [20-22]. This is still the 
case when the defect extends from a line in the bulk [23-27]. It was later shown, in 
a first order perturbation calculation, that the gap-exponent relation is valid for any 
marginal extended perturbation [28]. On the other hand, the geometrical shape of 
the free boundary may also lead to a modified critical behaviour. These effects are 
relevant in the critical behaviour at corners or parabola shaped systems (i.e. such that 
the boundary curve follows a parabolic law) [29-33] . 

In the present paper, we consider the two-dimensional Ising model with a marginal 
Hilhorst- van Leeuwen defect, as well as a related hyperbolic type of defect in the corner 
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geometry, and calculate the corresponding local critical exponents.. Using conformal 
methods, the problems are studied in the strip geometry. Here the spectrum of the 
transfer matrix is calculated exactly using the method of supersymmetric quantum 
mechanics. 

The setup of the paper is the following. In section 2, we present a short summary of 
SSQM and of the concept of shape invariance of the potential partners. In section 3, 
we show that in the cylinder geometry, when the Hamiltonian limit is considered, 
the eigenvalue equations in the continuum limit take the form of supersymmetric 
Schrodinger equations. The Hilhorst-van Leeuwen problem is considered and the 
complete spectrum of the transfer matrix is calculated by the method of SSQM. The 
same calculation is performed for the hyperbolic defect in section 4. In section 5, the 
critical exponents are calculated and a relation between the two problems through 
conformal invariance is discussed. 



2. Supersymmetric quantum mechanics 

The work of Witten [1] has focused considerable interest on supersymmetric quantum 
mechanics (for recent reviews see references [34,35]). Furthermore, by the concept of 
shape invariance, Gendenshtein [7] has obtained a systematic generalization to Dirac's 
operator method for the Id harmonic oscillator problem. 
Let us consider the Hamiltonian 

W- = -^ + v_(C) (i) 

with a vanishing ground state energy E . The ground state wave function is then 
related to the potential as V_(C) = V'o(C)/' i / , o(C)- I n terms of the superpotential 

W(C) = -^mVo(C), (2) 
the Hamiltonian 7iL is factorized: 

= + ( w2 (0 - = Q + Q~- (3) 

Here the prime denotes derivative with respect to ( and the charge operators are 
defined by: 

Q + = -^ + MC), Q- = ± + w(C). (4) 

The partner Hamiltonian 

H+ = ~ + V+(C) = ~ + (W 2 (C) + W'(0) (5) 
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may then be introduced and is also factorized, H + = Q~Q + , and there exists a one- 
to-one correspondence between the spectrum of the two partner Hamiltonians as: 
E~ +l — E+. If the ground state wave functions of H±, which are given by equation (2) 

as 



(0 = exp 



± J W(C)dC 



(6) 



are non normalizable, the ground state energies of both 7i_ and 7i + are non-zero and 
E~ = E+. In this case supersymmetry is broken. 

In the following we consider unbroken supersymmetry, i.e. Eq = and the 
potential partners which satisfy the shape invariance property as: 

V + ((,a )=V-((,a 1 ) + R(a 1 ). (7) 

Here a is a parameter of the Hamiltonian, a\ is some function of a , and R(a\) is 
a function which does not involve the variable (. It is then easy to show that the 
spectrum of 7i_ and H + are simply shifted by the amount of R(ai) and then, by 
iterating the shape invariance relation, one builds a hierarchy of Hamiltonians whose 
spectra are related as mentioned above. Finally, one finds the eigenvalues of H- as: 

n 

£-(a ) = 5>K). (8) 

k=l 

The corresponding wave functions are obtained by applying the charge operators on 
the ground state wave function: 

^n(Cao) ~Q + (a )(5 + (ai)---Q + (an_i)^o(C,an)- (9) 

The shape invariant potentials can be found in the literature [36-41]. The 
factorization technique was in fact originally introduced in the context of ordinary 
differential equations by Darboux [42-44], and the application of the so-called 
commutation formula to the Schrodinger equation can already be found in [45]. 



3. Hilhorst-van Leeuwen model 

Consider a semi-infinite two dimensional Ising model with inhomogeneous nearest 
neighbour couplings 

K(p,6) = K(oo)-gZ(p,0) (10) 

where K(oo) is the bulk critical value. The scale covariance requirement for the 
inhomogeneity leads to a power-law behaviour for the radial part of the shape function 
[28]: 



Z(p,9) = f(9)/p", 



(11) 
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and the perturbation amplitude g, then, scales under renormalization as g' = b yt ^g 
where y t is the bulk thermal exponent. Here, we use the method of conformal 
invariance. The deviation from the bulk coupling in the original system t(z) = 
K(p,9) — K(oo) transforms, under the conformal mapping w = w(z) = u + iv, 
according to t(w) =| w'(z) \~ yt t(z) [20]. With the usual plane-to-cylinder logarithmic 
conformal mapping w(z) = — In z, the semi-infinite system is mapped onto an infinitely 
long strip of width L with free boundary conditions, and the inhomogeneity (10) 
becomes 

K(u, v) = K'(oo) - g (£f exp [™(u, - y t )] f (^) , (12) 

where K'(oo) is the critical coupling in the modified geometry. If we furthermore 
assume a marginal inhomogeneity, i.e. such that the perturbation amplitude remains 
unchanged under a rescaling, one has ui = y t and it yields a perturbation which is 
independent of the u— direction along the strip: 

K{v) = K\oo) - 9 y (^-) . (13) 

The prototype of smoothly inhomogeneous systems has been introduced by 
Hilhorst and van Leeuwen [9] . Here, as an illustration, we recover the results previously 
obtained by Burkhardt and Igloi [20] by more complicated methods. Consider a two- 
dimensional semi-infinite Ising model on a square lattice. The couplings K\ parallel to 
the surface are constant, while the nearest neighbour couplings K 2 (y) perpendicular 
to the surface assume a power law deviation from their bulk critical value (figure (la)): 

K 2 (y) = K 2 (oo) - 9 -, (14) 

where y measures the distance from the free surface. This corresponds to a marginal 
shape function Z{p,9) = (psin^)" 1 . This model has been extensively studied in the 
two-dimensional classical version [9-16] as well as in its quantum counterpart [17,18,20- 
23] (for a review see reference [8]). Following Burkhardt and Igloi [20], we transform 
the inhomogeneity by the logarithmic conformal mapping and the inhomogeneity 
transforms into a sinusoidal form on the strip: 

K 2 (v)=K' 2 (oo)-y—^ 0<v<L. (15) 
Lsm(-) 

The transfer matrix along the strip T = e _TW , where r is the lattice spacing, leads, in 
the extreme anisotropic limit [46-48] , to a one-dimensional quantum chain defined by 
the Hamiltonian: 

L L-l 

n = °* W - o £ A(£)a^K(£ + 1), (16) 



2 e , 



with varying couplings 
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Here, the ex's are the Pauli matrices. The Hamiltonian H can be diagonalized 
by standard methods [49,50], transforming in terms of fermion creation (77^) and 
annihilation (rjk) operators as 

H = J2 A *Uvk-l). (18) 

Here, the fermionic modes with the lowest energies, which are 0(£ _1 ), are obtained 
in the continuum approximation from a pair of Schrodinger equations involving the 
inhomogeneity function x(0 = a / s ^ n C where ( = n£/L. The first one in terms of ipk 
reads as 

-^+(x 2 (0-x'(0)^(C)=(^)% fc (C), 0<C<^ (19a) 
with the boundary conditions 

^(C)\ -0 ^ 



= -X(vr). (196) 

C = 7T 



Similarly for the function <pk'- 



0UC) 



0*(C) 



+ (X 2 (C) + X'(O) 0*(C) = {^) MO, 0<(<n, (20a) 



= +x(0), fe (C)l c=7r = O. (206) 

C=o 



In these expressions, V'(C) an d 0(C) are the continuum limit approximations of 
the eigenvectors entering the discrete eigenvalue equations that one obtains when 
diagonalizing the Hamiltonian (16) (see reference [49]). 

The similarity between these equations and the Schrodinger equations encountered 
in supersymmetric quantum mechanics has been already mentioned by Choi [51], but 
here we show how the concept of shape invariance may be used to determine the 
excitation spectrum. 

First, we note that all the eigenvalues of equations (19a) and (20a) are the same, 
including the smallest one (A L/tt) 2 , thus, in the language of SSQM, supersymmetry is 
broken [52,53]. This statement is in agreement with Witten's argument [1], according 
to which unbroken supersymmetry requires a superpotential with one node (or an 
odd number of nodes). This is obviously not the case for x(0 ( see figure 2), 
which is symmetrical to tt/2, since the inhomogeneity in the semi-infinite plane is 
translationally invariant along the surface. We have then to face the problem of 
finding a superpotential W(£) in order to restore supersymmetry. This superpotential 
must be related to the inhomogeneity function x(0 by a Riccati equation, i.e. such 
that W 2 — W and x 2 — x' are identical up to a constant, the constant being of essential 
importance because its existence ensures that supersymmetry will be restored. 
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The boundary conditions in equations (19b) and (20b) generally pose the same 
requirement as in SSQM, i.e. the wavefunction must vanish at both ends of the interval 
since the potential term diverges there. However, if ijj' diverges faster than ijj when 
( — > 7r, then the solution of equation (19a) is non normalizable. This type of solution, 
which describes a localized mode, is associated to the appearence of spontaneous 
surface order in the system and corresponds to a vanishing excitation Ao = 0. For the 
Hilhorst-van Leeuwen inhomogeneity, such a solution is given by 



-0100.(0 ~ exp - / x(C)dC = tan 



(21) 



which is indeed non normalizable for a < —\. The lowest excitation energy in this 
region is then A = 0. We shall return later to determine the higher lying levels of the 
spectrum in this case. 

In the following, we deal with the region a > — |, where the method of SSQM 
works without limitations and shape invariance is a worthwhile concept to deduce 
the eigenvalue spectrum. First, we should find a convenient superpotential which 
solves the Riccati equation. This is done with the mapping introduced by Dutt et al 
[54]. The inhomogeneity function x(C) * s a special case of the Scarf superpotential 
S(() = a-i/ sin^ — a 2 cotC, which leads to the Eckart potential by S 2 (() — S'(Q [55]. 
With the choice a.\ = —1/2 and a 2 = a + 1/2, this defines a new superpotential 



w>(0 = - 



i 



2sinC 



a 



cotC 



(22) 



It is also easy to see that the superpotential W > (Q presents one node (figure 2), thus 
Witten's requirement on unbroken super symmetry is satisfied and super symmetry is 
now unbroken in the range a > — |. This choice leads to the trigonometric Eckart 
potential for V_(C): 



a 2 + acosC / 1 

V-(C) = '. 2~7 \ a+ o 

and the ground state excitation A can thus be identified as: 



(23) 



A n L 



7T 



a 



(24) 



Now equation (19a) can be written as: 



d 2 A 
dC 2 



+ (w 2 (0 -W>(0) MC) = 



A fc L 



7T 



A L 



n 



MO- 



(25) 



and the ground state wave function, obtained through equation (6) is given by: 

^(O-exp^-y W>(C)dc) ~-^sin a+1 C(l + cosCr 1/2 . (26) 
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The solution (26), which is indeed normalizable for a > —1/2, continuously evolves 
towards the localized mode (21) when a — > a c = — \ from above. The higher lying 
levels of the Schrodinger equation, which are given as 



AkL 



n 



AnL 



7T 



(27) 



are obtained from the shape invariance property of the partner potentials: 



V+(C,a ) = V_(C,ai)+ (ai + 



1 



a 



(28) 



where ao = a, a± = ao + 1. Then, according to equations (7) and (8), the 
energies of the single fermion excitations follow from the remainder function R(ai) = 
(a 1 + l/2) 2 -(a + l/2) 2 : 



7T 



A k = — [a + k+- 



k = 0,1,2 



a > --. 



(29) 



In the regime of surface order, a < —1/2, the eigenfunctions of the excited states 
of equation (19a) are normalizable, the previous method thus applies. Now, the 
superpotential is given by 



and the energy of the first non-vanishing excitation is identified as: 

2 /-, \ 2 



(30) 



AiL 



TV 



1 



— a 



(31) 



The higher lying excitations can be similarly obtained from the shape invariance 
property, so that the energies of the fermion modes are now given as: 



7T /, 1 
A = 0, A fc = - I k - a - - 



fe = l,2,3..., a< 



1 



(32) 



The potential V_ (C) and the corresponding eigenenergies E% are shown in figure 3 
in the ordered phase (a < a c ) and in the non-ordered phase (a > a c ) . Figure 4 shows 
the two first eigenfunctions in the two regimes. 



4. Hyperbolic defect 

The inhomogeneity in the Hilhorst-van Leeuwen model, studied before, can be 
considered as a result of elastic deformations on the free surface of the system. If 
the system has now the shape of a corner with a right angle, then uniform elastic 
deformations would result in a defect of hyperbolic form. The local couplings are 
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constant along the hyperbolas f(x,y) = = const, whereas the couplings pointing 
perpendicular to the f(x,y) lines are assumed to vary as (figure (lb)): 

K ± (x,y) = K ± (oo)-g-?-. (33) 

xy 

Thus, near to the surface but far from the corner, the inhomogeneity has the same 
shape as in the Hilhorst-van Leeuwen model in equation (10). The shape function 
corresponding to this defect is: 

ZM = T— ft - ( 34 ) 

p cos v sin u 

Once again, the inhomogeneity is mapped onto a strip geometry, now the appropriate 
conformal transformation is the the ^ In z logarithmic mapping. In the strip geometry, 
the inhomogeneity is again a sinusoidal form, the couplings in the Hamiltonian 
operator (16) vary as: 

W = (*t\ (*t\ =1 -I-hr (35) 

cos — sin — sin — 

\2LJ \2LJ \L 

from which we deduce the inhomogeneity function in the continuum limit: 

X(C) = ?. , = ~^r r (36) 
cos Q sin Q sin 2Q 

where ( is defined in the range < ( < n/2. This inhomogeneity function leads to 
the Poschl- Teller potential, but here, from the previous section we can immediately 
get the energies of the single particle excitations: 

A k = ^ z (2a + 2k + l), fc = 0,l,2..., a > ~, (37a) 
A = 0, A fc =^(2fe-l-2a), k= 1,2,3..., a < ~. (376) 



5. Local critical properties 

Conformal invariance makes it possible to transform critical systems from one 
restricted geometry into another, and deduce the local critical exponents in the former 
geometry from the energy gaps in the transformed one. 

In a semi-infinite system, like the Hilhorst-van Leeuwen model, the algebraic 
decay of the correlation function at the critical point between one point close to 
the surface {z\ ~ 1) and another point far in the bulk (z 2 ~ z) is asymptotically 
given as | z |~( x+x i) ) where x^ is the surface anomalous dimension of the operator 
fj,, while x is the corresponding exponent for the homogeneous bulk. Under the 
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logarithmic conformal mapping w(z) = ^flnz = u + iv, the correlation function 
transforms according to the usual position-dependent law involving only the bulk 
scaling dimensions: 

(fi( Wl )fi(w 2 )) =| w'( Zl ) H w\z 2 ) \- x (jiizM^)) (38) 

and the correlations in the cylinder geometry exhibit an exponential decay along the 
strip which defines the correlation length £ on the strip. In the extreme anisotropic 
limit, l/£ is given by the energy gap [47], so that the surface anomalous dimensions 
are contained in the spectrum of the Hamiltonian operator in equation (16): 

xt = -(E^-E,) (39) 

7T 

as L — > oo. Then, using the diagonal form of H in equation (18), the surface critical 
exponents of the Hilhorst-van Leeuwen model can be obtained as combinations of the 
A k fermion energies. For example the critical exponents of the surface magnetization 
and surface energy correlations are given by: 

1 1 

x? = a+-, xl=2a + 2, a > --, (40a) 



< = 0, x\ = \-a, a<~, (406) 

in agreement with [20]. Here, x™ = is due to surface ordering. 

For the hyperbolic defect, one defines the corner exponents, denoted x%, and 
associated to the algebraic decay of correlations in the corner geometry. In the strip 
geometry, the x^s are again proportional to the corresponding gaps of the Hamiltonian 
operator such as 

^ C = ^(E,-E ). (41) 

Then, the corner exponents for the magnetization and the energy for the hyperbolic 
defect with a right angle G = tt/2 are given as: 

x™ = 2a+l, x e c = 4a + 4, a> ~, (42a) 



x™ = 0, x e c = l-2a, «<-2' 

Comparing these results to those of the Hilhorst-van Leeuwen model, one can 
notice that the corner exponents are in each case the double of the corresponding 
surface ones. The same relation is known between exponents at a free surface and 
those of a corner of a right angle without the presence of an inhomogeneity, which is, 
according to Cardy [33] , a consequence of conformal invariance. The Schwarz mapping 
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z = z e ^ with = 7r/2 connects the two geometries and leads to the above relation 
between the local exponents. It is not difficult to see that the same Schwarz mapping 
transforms the Hilhorst-van Leeuwen inhomogeneity and the hyperbolic defect into 
each other, and thus gives the explanation for the observed relation between the 
corresponding local scaling dimensions. This last result can be used in the opposite 
direction, then the close relation between the spectrum of the Eckart and that of the 
Poschl- Teller potentials can be attributed to conformal symmetry. 

Finally we note that the relation between SSQM and inhomogeneous Ising models 
cannot be exploited further. Inspecting the table of shape invariant superpotentials 
[36-41], no further one is known at present which could serve as a basis for a new 
physically relevant inhomogeneity with an exact solution on the two-dimensional Ising 
model. 
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Figure captions 

Figure 1. Enhancement of local couplings near an extended defect, a) at a free surface 
(Hilhorst-van Leeuwen inhomogeneity) , b) at a corner (hyperbolic defect). 

Figure 2. Inhomogeneity function x(C) ( ) an d superpotential W(C) ( ) for the 

Hilhorst-van Leeuwen model for a = 2 (left) and a = — 4 (right). 

Figure 3. V_(C) potential ( ) and allowed eigenenergy levels ( ) for a = 2 

(left) and a = -4 (right) 

Figure 4. Ground state and first excited wave functions for the Hilhorst-van Leeuwen 
model for a = 2 (left) and a = — 4 (right). 



